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' The possibility that quantum corrections break the conservation of superhorizon adiabatic 

, perturbations in single field inflation is examined. I consider the lowest order corrections from 

' massless matter fields in the Hamiltonian formalism. Particular emphasis is therefore laid on 

I the renormalization. The counterterms are the same as in the Lagrangian formalism. The 

bX)' renormalized value of the tadpole is zero. I find a possible secular dependence of the power 

^ I spectrum at one loop due to the trace anomaly, but this result depends on the approximation 

■ of the modes and is inconclusive. The symmetry (not) violated by the quantum corrections is 

ly-^ I the invariance by dilatation. Perspectives on the backreaction problem are briefly discussed. 

(N \ 

Inflation would have probably remained at the stage of a beautiful theoretical idea if it had 
not been possible to relate cosmological perturbations between the inflationary period to the 
Q_^. matter and radiation spectra observed today. This possibility is truly amazing given that it 

D I implies relating cosmological perturbations between two periods separated by a vast energy 



range about which we know very little (reheating, phase transitions, and of course inflation 
itself). What saves us is the fact that the scales which are observed today where outside the 
horizon during most of the time that these phenomena took place, and were decoupled from 



> 

^ I what was happening on much shorter scales ^. 

'nI" ■ This decoupling takes on the form of a conservation law. It has long been known that 

, certain gauge-invariant combinations of the flrst order metric and matter perturbations become 

^ I constant outside the horizon This is the case for instance of the variable C = —ij) + -j^^^^^Siu 

■ where tp parametrizes (the diagonal part of the) perturbations of the spatial metric, H is the 

On . Hubble rate, po and po are the background energy density and pressure respectively, and 5u is 

^p. \ the potential of energy flow JT*^ = diSu. It is related to the intrinsic curvature of spatial surfaces 
by 7^ = —-^d'^C,. A mode of wave number q follows the equation 



where N is the number of efolds, i.e. a = exp(N). The flrst term (the non-adiabatic pressure 
Spnad = dp — ^5p) vanishes for superhorizon modes in many relevant cases, including single 
field inflation 0. The second term regroups gradients and vanishes exponentially fast once the 
mode exits the horizon. 

Two derivations of ([1]) are particularly interesting for their generality and the physical insight 
they bring. In the derivation offered in [1, 0], equation ([T]) and its nonlinear generalization 
descend directly from the conservation of energy projected on the hypersurface of constant 
time nP'VaT"'!, = (the vector n" is the unit normal to these hypersurfaces). In an alternative 
derivation Maldacena showed that in the infinite wavelength limit, the value of the action 
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By this I really mean that the mechanism at work is decoupling, and not acausality. By definition of inflation, 
the scale of causality is the size of the horizon at the beginning of inflation, and is therefore (much) larger than 
any observable scale today if infiation is to be of any use. 
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on a solution of the non-linear Hamilton constraint in single field inflation is a boundary term 

S = -2 [d^xdtdt (a^e^^H) + O f 41 (2) 



It does not contribute to the non-linear classical equations of motion, which therefore admit the 
solution ( = cte + 0{q^/a^). The physical contend of ([2]) and its relation to the first derivation 
of ([1]) is better revealed if we use the number of efolds = Ina as the time coordinate. One 
then notices that for q = 0, the action ([2]) changes by a total derivative when the coordinates 
are rescalled 

N + X, x' ^ e'^x' (3) 

Equation ([2]) thus implies the invariance by dilatation of the field equations at vanishing 
wavenumber, which is obviously a symmetry of ([T]). ^ 

So much for the classical theory. What can we expect in the quantum theory ? A breakdown 
of the conservation law would be tight to a violation of the invariance by dilatation. Now, 
under a Weyl transformation, i.e. a local scale transformation which leaves the coordinates 
fixed Qahix) I— > e^^^^^ gab{x), the effecive action obtained by integrating the matter degrees of 
freedom changes according to 

5nT = -\ [d'x^n{T^J (4) 
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The trace of massless fields is fixed by the curvature, (T*^) = 0{H^). (This is true whether 
the field is conformaly coupled or not, and by abuse of language I will call this term the trace 
anomaly regardless of the coupling to gravity.) We can conclude that if non-conservation of C 
there is, it is because of the trace anomaly. This is only a heuristic argument. Whether this 
violation actually happens depends on the details of the interactions (because the operator Tab 
is inserted into loops and is therefore integrated over time with a certain weight function), and 
on the details of the renormalization. This is the calculation reported in this paper. 

The rest of the paper is more technical. The next section motivates this calculation and 
presents its technical aspects which are expanded in considerable details in the subsequent five 
sections. The reader who does not want to get involved in the technicalities is invited to read 
the next section to get an idea of it and proceed to section IVIII where the main results are 
summarized and discussed. 



I. INTRODUCTION TO THE TECHNICAL ASPECTS OF THE CALCULATION 

The conservation law is expressed as the time independence of the power spectrum of super- 
horizon modes: 

V{q,t)^ ld'xe-'^^{C{t,^)at,0))^V{q), for (5) 

A first important step was achieved by Weinberg 0, 0]. Using the Hamiltonian formalism 
he showed that in a large class of theories, the late time contribution of quantum loops ^ is 



^ The fact that a rescahng of the conformal time corresponds to a translation of the variable A'^ makes it intuitively 
clear that the two derivations are related, since invariance by translation is associated with conservation of 
energy. To be more precise, there exists an equivalence class of coordinate systems in which the metric reads 
ds^ = a^e^'' {—dr^ + dx'^) up to corrections 0{q^ /a?) This implies that the normal to the constant time 

hypersurfaces n becomes aligned with the conformal Killing vector d/dir. 

Unfortunately cosmologists use the word 'loop' to describe the perturbative corrections in the so-called 5N 
formalism. These corrections mix the quantum mechanical correction to a purely classical one which accounts 
for the passage from the flat gauge to the comoving gauge. These are not the loops considered here. 
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at most powers of ln(a). The proof of this theorem, which is also the way calculations in 
the Hamiltonian formalism have been done ever since, rests on the possibility to exchange the 
order of integrations, integrating first over time at fixed external wavenumber q and internal 
wavenumber p, and to take the limit t — > oo before integrating over momentum. The integral 
converges at early times because of the oscillatory behaviour ~ e"*^"^ of the mode functions, and 
in the conditions of the theorem it also converges at late times because the integrand contains 
sufficient factors of l/o. The remaing integrals over the loop momenta are easely evaluated in 
dimensional regularization. This exchange is of course only allowed for convergent integrals, 
which means that the ultraviolet divergences must be cancelled by counterterms. If such is the 
case, the integral over momentum is effectively dominated by values p ^ q, and the value of the 
whole loop is effectively set at ~ — 1 like tree graphs. 

An important point is that the ln(a) do not necessarily occur. The proof of the theorem is 
algebraic (one counts the number of factors of a), and the theorem only states that they can 
occur (as the limiting case). Whether they do or not depends on the details of the interactions 
and on the renormalization. To clarify these subtleties is one of the objectives of this work. 

Why should one care about this, and why is it a non trivial problem ? The first question is a 
legitimate one indeed. The possible violation of the conservation of C is certainly not an urgent 
problem because if it occurs, it is likely to be too small to have observable consequences. My 
answer to the first question is that the calculation reported here should be considered as a baby 
step towards a better understanding of the backreaction problem, in two respects: a theoretical 
one, with regard to the identification of proper observables to quantify the effect, and a pratical 
one, to develop and test the techniques to be used in future calculations of the backreaction. (I 
will return to this in section IVlIFi ) Indeed, the Hamiltonian formalism lends itself particularly 
to calculations of the correlation functions of cosmological perturbations because one can work 
with the physically relevant (non linear) variables, e.g. C ^-nd the two helicities of the graviton. 
Follows a considerable simplification of the algebra compared to the Lagrangian formalism where 
calculations are carried out with a 4 x 4 tensor hat (the 'pseudo-graviton') and involve ghost 
fields. The downside of adopting the Hamiltonian formalism is that by abandoning manifest 
Lorentz invariance, we are deprived of a powerful instrument to identify counterterms. This is 
the answer to the second question. But since we have very good reasons to use the Hamiltonian 
formalism, we must try and make it work. 

The task consists of the following steps which also delineate the plan of the paper. First, 
the Hamiltonian is calculated up to fourth order. I will ignore gravitons and I will only calcu- 
late the corrections from matter at one loop order. These assumptions considerably simplify 
the calculations which are already relatively tedious. The part of the interaction Hamiltonian 
relevant for this calculation consists of two pieces, a three vertex ~ (^aa and a four vertex 
H4 ~ CC^f^- The calculation of H4, which is the first new result of the paper, is a bit lengthy 
because the constraints must be solved at second order and the relation between field derivatives 
and momenta is non trivial by virtue of the derivative couplings between matter and gravity. 
The canonical quantization is detailed in sec. Ill D[ 

Before calculating the loop corrections I present in section IIIII the expressions of the coun- 
terterms quadratic in the curvature, namely the Weyl tensor = C'^^'^'^Cabcdi the Euler density 
E4, the square of the Ricci scalar and the Dalembertian of the Ricci scalar. 

Section HV] is devoted to the renormalization of the tadpole {0, which is tantamount to the 
renormalization of the energy-momentum tensor. It is shown that the renormalized tadpole can 
be consistently set to zero. 

The corrections to the power spectrum from H3 and i/4 are considered in sections |V] and sec. 
IVII They are shown to be renormalized by the same geometric counterterms as the tadpole. 
Several appendices complete the presentation of these calculations. The following five sections 
are very factual. The important results are discussed in sec. IVIII 
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II. CANONICAL QUANTIZATION 
A. Model and conventions 

I use the model of 0]. Matter consists of one inflaton field (p{t,:x.) characterized by a classical 
potential V{ip) and a vector field (?(t,x) = "^(cJi, ...) in the defining representation of 0{J\f), 
with no self- interactions. ip{t,x) and (j{t,x.) are indirectly coupled through gravity. We work in 
the gauge where the inflaton is homogeneous on each space-like hypersurface, 

(^(t,x) =(^(t), (6) 

and where the parametrization of the spatial sections exhibits the physical degrees of freedom 
of the problem 

hij = a^{t)e^'^ (exp7)jj , 7m = , daij = . (7) 
where a{t) is the background scale factor. The remaining metric elements are 

900 = -iV' + h,jN'W , goi = h,,N^ (8) 

where and A^* are the lapse and shift functions. The comoving time is noted t = J dt/a{t). 
The Lagrangian density is 

C = (i? + - 2V{ip) + g^'daadba 
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H — |^(T — N^di(jj — [exp (—7)]*-' diadjfj + total derivative (9) 



The second form is the one appropriate for the canonical formulation [3]. The new terms 
appearing in this equation are: 

Eij = ^ (dthij - V^Nj - VjN,) (10) 

related to the extrinsic curvature of constant time hypersurfaces by Kij = N~^Eij, its trace 

E = W^Eij, and the Ricci scalar curvature of spacelike surfaces TZ = W^TZij. I will ignore 

gravitons. Since the cjj are essentially N copies of the same field, I will be sloppy in my notations 

and often write a for a. 

Before proceeding with the solution of the constraints and the canonical quantization, I 

introduce a few notations and conventions. I will from now on write the scale factor a{t) = 

■ 2 

and the Hubble rate p. The first slow-roll parameter is defined by e = — ^ = To simplify the 
Fourier transform of the upcoming expressions, the spatial gradient di = and Laplacian are 
defined with respect to the flat metric 5ij. Hence 9* = dj = di, and = q^t^ = ^^•'didj. 

The inverse of the Laplacian is noted d~^. I will ofen use the notation d^p for the measure ^^^^3 . 
Conventions for the curvature are those of the book of Misner, Thorne and Wheeler The 
components of the Riemann tensor are R^'f,^ = dcT'^^ — ddTf^ + Tg^F^^ — F^^F^^, and the Ricci 
tensor is Rm = R%ad- 

The algebra will be done in full generality, but for the calculation of the integrals in closed 
form some approximate solution of the mode equations must be chosen. The linear equation of 
C is (see the quadratic part of the action 

C+Up + -]c-^d^C = o (11) 
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The linear equation of a massless scalar field with coupling i^—gRa"^ is 



2 

a + ?,pd + \2i(?a - i- = (12) 

Following I will use the quasi-de Sitter approximation around the time of horizon crossing. 
The positive frequency solutions in the asymptotic past (i.e. corresponding to the free Bunch 
Davies vacuum) are then 



C,(t) = C° (1 + ^qr) e-"- , \cf = ^ ^ (13) 
For a minimally coupled scalar field = we have 



aq{t) = (1 + iqr) e^^^ , \af = ^ (14) 



and in the case of a conformal coupling ^ = 1/6, 



(J° 1 



-, = f^-. KT = ^ (15) 

The normalization constants C,^ and are fixed by the canonical commutation relations. 

I will discuss in due course the validity of this approximation, namely in sections WN\ IVII El 
and appendix O 



B. Solutions of the energy and momentum constraints 

The Hamiltonian and momentum constraints obtained by variation of ([9]) w.r.t. and N"^ 
are respectively 

iV^ {TZ-2V - h'^diadja) = EijE'^ - + ip^ + (^^ - N'dia^^ (16) 

[N-^ {E) - 6'jE)] = ^djaia - N'd^a) (17) 

They are solved iteratively using the Ansatze 

N = l + a + a2 + N' = d,{B + B2)+w' + wi + ..., diw\^) = (18) 

I will write first order quantities without an index to avoid an exponential number of occurrences 
of the number one. Some expressions are simpler if we use the conformal time r and I will go 
freely from t to r. For the latter, it proves convenient to absorb the factor of a from dt = adr 
into new parameters 

uj = aB , = aw'' (19) 

so that we have hijN^dtdx^ = a? e~'^'~- {diB + w^)dtdx'' = e~'^^ {diUJ + w'^)dTdx^ . 

Only the first two orders are necessary to calculate the action up to fourth order H (see 
footnote H]) . The solutions at first order are 

a = - = ^ (20a) 
P P 

B = -^ + ed-^C, LJ = -^ + ed-\' (20b) 

w' = (20c) 
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d'^am = -^^djiadja) (22) 
4p V «W P 



They do not depend on a because the field has a vanishing v.e.v., and they do not depend on 
the inflaton fluctuations because we work in the foliation Q. At second order, the lapse and 
shift receive separate contributions from matter and first order metric perturbations, 

"2 = + , B2 = Bg + Bm, W2 = w'g + (21) 

The purely matter contributions are 

1 

d'Bm = -^[^'' + ) - -"m (23) 

d'^wi, = -2adja + ApdjUm (24) 
and the contributions from first order metric perturbations are 

2pd'^ag = -292 (diCdiB) + dj { {djad'^B - diadidjB) + SdiCdidjB - djCd'^B} (25) 
-^d^wi = -2^J[p^,ag + {^,C^^B)] 

+ {djad'^B - diadidjB) + SdiCdidjB - djCd^B (26) 

and 

Ap'd'^LOg = -Aa'^Vug - 6eC'^ - 2{dCf + SCd'^C - l2p'diCdiUj + [{d'^ujf - didjtodidjuj] (27) 

C. The action up to fourth order 

We can now expand the action ([9]) up to fourth order ^. The quadartic part is 

s,.iL^!,jc^-mu^^-m] (28) 



2.1 IV a2 

The cubic interactions are 

\2 /3 



S3 = - [aU -2eC-^^ - 2e^ - Ad^BdiBdiC - - W^B did^B - (d'^B) 
J [ « P p ^ ' ' 



I substituted C/p ^or a but left the other parameters of the lapse and shift to simplify the 
expressions. The first line of 5*3 regroups the purely gravitational cubic interactions CCC (since 
B is independent of a) and the second line regroups gravity-matter interactions (aa. This result 



"* I brielfy recall why the solutions of the constraints at third and fourth order are not needed for this calculation 
[^. Upon expanding the action ([9]), one observes that 04 multiplies the zeroth order energy constraint, and 
as multiplies the first order energy constraint, which are both assumed to hold; terms in VA'^ of third and 
fourth order cancel; V(iA'j) V'-'A^-'^ starts at second order, but after integration by parts the second order term 
turns out to vanish. The final form (|28I29|I is obtained after several integration by parts, use of the constraints 
equations, and having dropped spatial divergences. 
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is not new: (^^ interactions were derived in 0] and the matter-(^ interactions in 0|. The fourth 
order interactions are 



+ 



diw{ diw{ - Ad'^B diB2 diC - Ad'^B2 diB diC - 6{diB diCf 



2 '^t'^2 

i 



2didjB didjB2 - 2d'^B d'^B2 + 2didjB diw{ - Ad'^B diB diC 



- 2&diB2 dia - 2&wi d,a + {^^B diuf - + Cjj 



The pure C part was previously given in ^ 

Only the Hamiltonian of the (^aa interactions has been previously calculated 0. I present 
the derivation of the complete Hamiltonian at fourth order in sec. IIIDI Before this, I check the 
expression of the matter-i^ interactions against the action given by the covariant formalism. In 
terms of the pseudo graviton tensor hat defined by 

ds'^ = (r) {T]ah + Kh) d.x"-dx^ (29) 

the action up to fourth order is 

S=-\ Jd^x^g'^'daadta = Jd'xa^r) (h^' + H^'^ dacrd.a (30) 
where the two auxilliary tensors are 

Kb = Kb - ^Vab , h = rf'-Kd 

H-" = Kch\ - \hKb + ^ (x - ^-'^^"') (31) 

Their indices are raised and lowered with rjab- Substitution of ([7]) and the solutions of the 
constrains into (f3T]) gives 



Sr = \ja'\^i^C-y^a^- 2d.Bad.a - + ^ ^| (32) 
which is equal to the second line of (f29l) . and 

Sr - \ ja^ [-a2 (.^ + ^) + + le - 3C^) - 2 (3C - ^) aOM. 

2 (e , .C^ {daf^ 



2&diB2dia - 2&w'dia + {diBdiaY - V + 



2 p j a 



2 



For comparison with this reference, note the use of different conventions: — e ^^Sijdidj and (9C)^ = 
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But for the first term proportional to 02, this is the last two lines of (|29p . There is actually no 
discrepancy between the two results: the terms — Q2 ^o"^ + ^^^] ^ in (j33p can be combined with 
the other terms proportional to 02 and simplified using the second order constraint equation to 
yield the term 2V{a2)'^ that appears in ([29|) . 



D. Canonical quantization 



I now proceed with the canonical quantization. One first calculates the conjugate momenta 

5S 

., = ^, vr. = - (33) 
and inverts these expressions. Schematically we have 

C = C(vrc, C, T^a,(T) =TT^ + 62C + S3C + ••. (34) 

a = fT(7r^,C,7r^,o-) = tTct + (52fT + (53(7 + ... (35) 

where 62 (resp. 6^) regroups the terms second (resp. third) order in the field variables, and the 
dots represent higher order terms. To simplify the expressions, I do not write the factors of a'^ 
and e associated with each factor of vr. One then substitutes these expressions into 



H = j(fx {vr^C + -^a^ - {C,2 + £3 + A + •••)} 



(36) 



In the fourth order Lagrangian £4 we can directly replace Q by ttq and (T by vTo-. The substitution 
in £3 requires the expressions ()34p up to second order in the field variables, and we a priori 
need these relations up to third order to substitute them into vr^^C + '^o-o' — £2- It turns out that 
these contributions cancel: 



^cC - 'C2(C, C) = (^c + -^aC + 
_ 1 

~ 2 



- {S2Cf + 



(37) 



and similarly for the (j-sector. The dots stand for higher order terms. Thus we only need to 
calculate the momenta at second order, a most welcome simplification. 

Because of the non-local form of the interactions, it is easier to proceed with the Fourier 
representation of the action. With the convention that repeated indices are summed, the cubic 
action has the form 

S3 = ^5(p + k + l)|cp,k,iCpCkCi + rfp,k,iCpCkCi + ep,k,iCpCkCi} 

r/ 1 



+ -5(p + k + l) 



3Cp 



P 



0"kO"l + /p,k,lO"pCrkCl + fi'p,k,lO"pO"kCl + V.k.lCpO'kO"! { + ■■■ (38) 



The dots stand for the terms without time derivatives. The expressions of c,d, ... are read off 



Cp,k,l 
(^p,k,l 

ep,k,i 

/p,k,l 



- I 2 + e 

P 



-2e 



<T(k,l) 
k2 12 
1 <T(k,l) 

a^p k2 
cT(k,l) 



+ 2e- 



kl 



k2 



2e d 



p,k,l 



a2/j (ap)2 ~'~ p2^ 
kl kl 



a2/9 



5p,k,l 



-2e- 



12 



¥,k,l 



k l 

a2p 



(39) 
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and the auxiliary function is 

a{k,l) = {k-lf -kH"^ (40) 

Note that these vertices have not been symmetrized, so that it is important to respect the 
ordering of the momenta according to the definitions (j38p of the functions c,d, .... 

We wiU later need the value of these coefficients for p = 0, k = —1. Most of the terms 
above are singular because they originate from the inversion of the solution d^B = and this 
equation is not defined for p = 0. This quandary is resolved once one remembers that the zero 
mode of C is a gauge mode and as such does not posses a well defined momentum. So we must 
for consistency take the value of these coefficients to be zero. 

Taking the variations of ()38p with respect to Cq and (Tq and inverting the expressions up to 
second order in the field variables, one gets 

hcr-c^ = - ^<5(q + p + k) 1 2 ^3Cp - vr^^ + /q,p,k<7pCk + 5q,p,kfTp7rli^| (41) 

2e (52C-~q = -^'5(q + P + k) |Cq,p,k7rlp7rij^i)q,p,kvr£pCk + eq,p,kCpCk} 

-^'5(q + P + k) |-^vr!;p7r:ik + 5p,k,q7r1pCrk + /iq,p,kO-pO-k| (42) 

where 

C'q,p,k — Cq p k ~\~ Cp q k ~l~ Cp k,q 

-Cq,p,k = f^q.p.k + ^^p,q,k (43) 

I recall that to simplify the expressions, I have not written the factors of and e on the r.h.s. 
The expression of 52Ct has a simple form in the position representation 

S2&ix) = - ^3C - ^ + N'dia (44) 

which could have been obtained directly from Q. The expression of 52C{x) is on the other hand 
far less neat-looking and I do not write it. Substitution of these expressions into (|36p and using 
the expressions (|4m42p to simplify some terms yields 



H3 = - i'lCi (46) 



2e{62Cy + {hay (47) 



The + sign in front of the last two terms of (jSOp is not a typo. It is the sum of —^(^20^ 
from (j37p and +((52C)^ obtained by collecting the terms produced by the variation of £3 in the 
substitution C *■ tt + 62C (by the very definition of 52C)- These last two terms are produced 
by the inversion of the relation 7r(^) and are therefore peculiar to the canonical quantization of 
theories with derivative couplings. 

The procedure of canonical quantization is completed after replacing ((^, vr) in the above 
expressions by the interacting picture variables (C/, vr/). Their evolution is by definition governed 
by H2, and the interaction Hamiltonian is given by (j46p and (|47p with the replacement 

7r<^ I— > a^e C/ , vTo- 1— > a^(T/ (48) 
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In the following I do not write the indices I. I only write the part of the interacting Hamiltonian 
relevant for the calculation of matter loops: 



a 

Y 



3C 



c 



(49) 



and 



e , X\ {da) 



(50) 



The first vertex is the sum of the relevant part of —£4 and ^{52^)"^. The second term stems 
from having fixed the gauge and solved the constrains 



+ ( 3C - ^ ) [2didjB didjBm - 2d^B d'^B^ + 2didjB d, 



(51) 



The third term in (j50p comes from the inversion of the relation between time derivatives of the 
C and ttq, i.e. the last two terms in (j47p . In the momentum representation, it is given by 



rr extra 
^4 



a 



X (5(p + k - 1 - I' 



11' r • q . 

— i^cTiay — 2e — :^d-\a\i 



1 



Cq,p,kCpCk + ^q,p,kCpCk + ^ (*^q,P.k + e-q,p,k) CpCk 



(52) 



These additional interactions are the price to pay for the tremendous reduction of the number 
of variables compared to the covariant formalism. One way or the other, appearance of such 
terms is inevitable because we need to fix a gauge in order to quantize the theory, which in the 
covariant formalism comes at the price of ghost fields. 

The basic observation one can make at this point is that the passage from (|33p to (|50p 
generated a certain number of new terms. It is thus not obvious that the ensuing divergences 
can be absorbed by local geometrical counterterms as in the Lagrangian formalism. It turns out 
however that both H^'^ and Hf^^^^ have a vanishing contribution at one loop. 



III. GEOMETRIC COUNTERTERMS 

When I carried out the canonical quantization I did not include the counterterms in the action. 
I justify this choice by referingto the point of view of effective theories, which implies a consistent 



perturbative treatment 10|, lul] . Because the effective action is a perturbative expansion in /i, 
counterterms of dimension [mass]" with n > 4 produced order by order by the loop expansion 
must not be considered as part of the classical action, but only as additional vertices. This means 
in particular that they must be discarded from the propagator (in the Lagrangian formalism) 
and from the canonical momenta (in the Hamiltonian formalism), despite that they contain 
terms linear and quadratic in C,. I will therefore take j dt Hcnt = — •S'cnt- 

The expressions of the Weyl tensor to first order and of the Ricci tensor and Ricci scalar to 
second order can be found in appendix [Bl 
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A. Squared Weyl tensor 

Let us begin with the squared Weyl tensor, 



= ^ jd^'x^C^ = ^ jdTd^x32{ep'Cf + 0{e) (53) 
As I explained I posit that the counterterm in the Hamiltonian formalism is 

= -mC^iep'f I ^ C;Cp (54) 
With the notation Q{t) = ^q(T)C-q(T"), the corresponding correction to the power spectrum is 

AwV = -2Im j dn (Fw(Ti)Q(r)) 

J —oo 

= 64Cwe^/m| T dr, [p'C^(ti)C*(t)] ' 



= 32Cw(6|C,T) 9V) (55) 

In sec. |V] this term will be shown to absorbe the divergences of the corrections from the 
interaction C,a(j. 



B. Squared Ricci scalar 

The other counterterms are more involved because the background is not Ricci flat, so they 
must be calculated at second order. The simplest is to start with the Ricci scalar squared. Let 
me note ^ = 6^^ = 6/)'^ (2 - e). We have, 

Sr^ = ^ jd'x^gR' (56) 
= ^ y"i?2 j^l + 3C + « + ^C' + 3C«) + SRep'C + r(^{2 + e)T,, + ^ 
+ AR {^^{d'^C? - ad'^a + (1 + 5e + £3) ad'^C - (3 + 2e) C^^C^ 
+ 432p'\c' + 624e(/>'C')^ - 296(ep'C')^ 

Here, Tqo = (o"')^ + (dcr)'^- I neglected terms second order in the slow roll parameters e = 1 — j^ryi 

and €2 = y^- Let me comment the various terms appearing in this expression: 

1) The list of quadratic terms is exhaustive: no other term quadratic in ^ with four derivatives, 
including a maximum of two time derivatives of C, can be constructed. Therefore, the other 
counterterms R"'^'^'^ Rated, R%I^a ^'^d Di? contain the same list of terms, only the prefactors 
differ. This is checked explicitly in appendix [Bj 

2) One of the quadratic terms in (|56p is (ep'C)^, like C^. There is therefore a degeneracy between 
those two counterterms that cannot be lifted from the sole term {HaHaCO of sec. |Vl 

3) Appear also terms quadratic in the matter fields. They originate from the solution of the 
shift vector ([23l) . d'^LVm = + .... They give a vanishing correction at lowest order 

2Im J dn (FcntCq(r)C-q(r)) = 12|C,(T)|2/,n J dn ^(ri)(-4(a')' + (2 + e)roo) = (57) 
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4) The quadratic terms (|C^ + 3(^0;) come from the expansion of \/—g. They give a vanishing 
correction. 

5) The terms p''^C,C,' give a correction 

APoc-2/m fdnp''\CiC'iC\r)) = V | (cOC*(r))' (1 + zgr^e-^gn | 

= - V|qO|4 {In(-gT) + cte + 0(gr)} (58) 
I used the asymptotic behaviour of the exponential integral 

^-igri ^ in(-gr) + C + 0{qT) (59) 



It is the only quadratic term of (j56p with this logarithmic dependence. It will find its use in 
section IVI CI The other quadratic terms of ()B14p containing (^'^ or gradients give the same 
polynomial corrections: 

per, I^Cd\, l^'^d'C, /^(^'C)^ l^ad'a ^qX\' (60) 

These terms will be used to renormalize the corrections of the fourth order interactions. 
6) Finally there are terms linear in ^. They appeared because the background is not Ricci flat. 
They serve to renormalize the tadpole As we will need other counterterms in addition to 
the Ricci scalar squared, let me calculate the contribution of a general linear counterterm 



Hun = -^, jd'x (37C + 5^) 



(61) 



A ^ 


^^ = 3^' 
6 . 

TGat = 


<5a=2 


Gn 


Sgn = 




7R2 = 36(2 - ef , 


Sr2 = 


E4 ^ 


7E4 = , 


<5e4 = - 


DR 


7n = 0, 


<5n = 4 



where 7 and 6 are constants. This expression also includes the linear terms from counterterms 
parametrized by a cosmological constant A f \J—g (or more generally of the inflaton potential) 
and a shift of the Newton constant J yJ—gR. For these terms as well as for the squared Ricci 
scalar, the Euler density E4 = R"'^'"^Rahcd — 4i?"''i?ab + -R^j and Di?, their values are 

1 

4(6 -e) 
16e 

(62) 

The corresponding correction to the tadpole is 

A(Cq(r))cnt = -2Im!^J dri{HnnCci{r)) 

= 5(q) X (5 - 7) qX\^ { Ini-qr) + cte + O^qr)} (63) 

Note that although neither 7^2 nor (5^2 are of order e, their difference is. This will prove 
important in the following section. Similarly for G^. We also see that the cosmological constant 
has a vanishing contribution. 

In conclusion, it is a significant check that the counterterms are of the same order in e as 
the interactions they renormalize: H3 induces corrections of 0(e) to the tadpole, which are 
renormalized by the linear terms ()6ip . and corrections O(e^) to the two point function, which 
are renormalized by the square of the Weyl tensor, while H4 induces corrections 0(e'^) which 
are renormalized by the counterterms ()58p and (I60p . 
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IV. RENORMALIZATION OF THE TADPOLE AND TRACE ANOMALY 

The renormalization of the tadpole is tantamount to the renormahzation of the energy mo- 
mentum tensor since = jd^xh°'^Tah- Indeed we have 

(Cq(r)) = 5(q) X 2Im |c*(r) dn a\ (-ai(Too) + Ci5*^ (7^,)) | (64) 

where 2Too = cr'^ + (9(t)^. We immediately see that the approximation is inadequate. If 
one insists on using it, one finds after integrating over time 

(Cq(r)) = 5(q) X {^(1 - er^)-h - (65) 

2 

where Jq and J2 the g-independent divergent integrals such that Too = Jo + "y<^2- There is no 
term of the list (|56p with the proper g-dependence to absorbe these divergences. 

Of course we know how to proceed. The renormalized value of the energy- momentum tensor of 
fields in curved backgrounds is obtained from the adiabatic expansion of the modes (the covariant 
methods look very different but the requirement of the Hadamard form expresses essentially 
the same idea. I therefore favour the more intuitive language of the adiabatic substraction). 
Since this expansion is close to an expansion in slow roll parameters, we understand that the 
instantaneous de Sitter approximation around horizon exit (jl4J) cannot do the job. 

The counterterms for a minimally coupled scalar field are [191 ] 

^rcnt = -f^^T^ {3C2-E4-12ni? + 5i?2} (66) 

(The last two terms are absent for a conformally coupled scalar field.) I make a simplification 
and take the renormalized value of the energy momentum tensor in de Sitter space 

{Tal,) = id gah (67) 

where /3 7^ is a constant. The trace (T"^)i.en = 4/3/9^ therefore does not vanish, even for 
conformal fields. For minimaly coupled fields in the Allen-Follaci vacuum, /3 = After 
substituting (f67l) and integrating over time we get 

(Cq(r)) = -/?<5(q) x/m|c,*(r) (3Ci+ai)} 

= (68) 

This is a significant result because the symmetry responsible for the conservation of ^ is a 
symmetry of the action for q = 0. True, the zero mode of C is a gauge mode. But like a 
Goldstone mode it acquires a physical relevance by the fact that it can be continuously extented 
to physical modes [13]. Had we any other linear combination than 3C + a in ([68]), we would have 
obtained (Cq=o(''")) oc ln(— g'r) (as in ([63]) ). This perfect fit of coefficients bears no relation with 
a fine tunning. It is a sign of the consistency of the theory, to wit, the symmetry fixing (1670 is 
not responsible of its own demise, as would constitute a breaking of its subgroup ([3]). 

Two other results back up this assertion. One is that the corrections from the linear coun- 
terterms (|62ti63p are 0(e), and therefore vanish in the de Sitter limit. This is consistent with 
the fact that we are neglecting interactions and that linear scalar perturbations in de Sitter 
space are gauge modes. The second result is that if we do not assume the de Sitter symmetry 
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this cancellation does not happen. In a model of power law inflation a(r) = r ^ " with a > 
one obtains for qr < 1 (see appendix [C]) 



(Cq(T)) = #<^(q) X q^qrY 



(69) 



which is finite for a finite value of r and — > in the asymptotic future. Moreover, the eln(— gr) 
in (|63|) is also replaced by {qr)", which shows that the linear terms (|61|) are of the appropriate 
form to renormalize the tadpole. 



V. CORRECTION FROM TWO INSERTIONS OF H3 



Instead of calculating the corrections from (09]), it is possible to make a cunning canonical 
transformation such that the scale dependent correction q^ln{q/fi) has a unique antecedent 0| 
(the corresponding transformation for the interaction is given in [^). I quote the result. 
Using again the notation Q{t) = (^q(t)C-q(t), we have 



AV 



t rt2 

dt2 / dh { Hsih), Hs{t2), Q{t) ) 



dtl { 



, Q{t) 



The new vertices //•^ and are 



Hs = Jd'x (Too + n) {-epa^d-\ 
T = cC' fd^xBToo 



(70) 

(71) 
(72) 



One clearly sees from ()70p that only (j7ip can be responsible for a q'^ln{q/ij,) because the term 
(j72p generates effective fourth order interactions which produce polynomial corrections g^. 

I now make a remark that will shortly be useful: one recognizes H3 as — ^ fd^x^/—g 5gabT°'^ 
in the longitudinal gauge, where the gravitational potential (pi and the perturbation of the three 
metric ijji are equal and given by 



(1)1 = ^^1 = -epa'd-^C 



(73) 



The proof of this statement can be found in appendix Rl 

I now proceed with the corrections from (j7ip . The finite part has already been calculated 
in 0], so I focus on the renormalization. To show that the divergence can be absobed in the 
counterterm , I exploit the underlying analytic and tensor structure made obvious by the fact 
that the Hamiltonian is — Jd^xC^. 

After doing the Wick contractions, the first term of (j70p can be written 



X / dT2V2 r dnViReiz^q'^Mlp, 
J —00 J —00 



(74) 



is the number of scalar fields. The vertex function is V{t) = —eHa^. The function Zq comes 
from the Wick contractions of (CiC2Ct^) — {CiCtC2)- It is given by 



Z 



1 



q'^aia2 



a(ri)C*(r) [C',iT2)C:{r)-UrK{T2)] 



(75) 



15 



7W comes from the Wick contractions of ((T^(ti)iT^(t2)), 

2 



-^PP'(n,r2) = -^a'^{n)a';{T2)a'p,{n)a'/{T2) (76) 
If in (j74p we factor out the normahzation of the metric modes we obtain 

(Cq(t)C-q(t)) = -IGATellCX [ H Re {(1 - iqrf S> - (l + gV^) S<} (77) 



where 



E> = e-^(-^+-"2-) 1^ (2^)^5(3) (q + p + p') A^,ed(p,p', n, r^) 



(78) 



S< = e-^(--^^) y^|!^(2vr)35(3) (q + p + p') Alrcd(p,p', n, r^) (79) 



and 



r, g-i(p+p')(Tl--r2) 

M,,d{p,p',n,T2) = {pp'Y (80) 

The latter is indentical to the positive Wightman function ((^^^^(''"i))^ (f5r2X('^2))^) of a scalar 
field X in Minkowski space in its ground state. We notice that the factors of a and p cancel. 
It is a simple exercice to verify that the expression (j74p is the Fourier transform of 

AGF(r,r',x - y) = iM [ dT2(fzi f dTid^Z2 (G>(r, y; r2, Z2) - G<(r, y; r2, 



00 J —00 



Gl.(T,x;Ti,zi)i;^(ri,zi;T2,Z2) - G^(r, x; ri, zi)I;^(ti, zi; r2, Z2) (81) 



where Gi^ is 



Gi(r,n) = Gl(n,r) = /(^a e-(-)C,(t)^ 

= /(03^^'^^^-^^K,T(l + ^9r)e-''(--) (82) 

It is the positive Wightman function of Q with the substitution Cg(^i) Cq{ti)/Q^ foi" the leg 
going into the loop. The self-energies are the connected part of the fluctuation of the energy- 
momentum tensor of the scalar field 

S^(ri,r2) = S^(T2,ri) 

= W{n)V{T2){a\h)a''{t2))cor. 

= Ae\lp,alp2 /|^e^(P+P')(^-^^)-Mred(p,y,ri,r2) (83) 

where A^red is defined at eq. ([80j) . The conjugation relations G>(r, ri) = [G<(r, ri)]* and 
similarly for show that the r.h.s. of (jSip is real. 

This game of rewriting has a purpose. We have put eq. (j74p into the generic form (|8ip of 
the one-loop correction of the two-point function in the in-in formalism (see Appendix [D] for 
an outline). Moreover, we see on (j80p . ()82p and ()83p that the various factors of o and p cancel 
before integration and the only time dependence is through oscillating functions. This is an 
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important remark since it implies that, as far as the analytic structure is concerned, except for 
the external legs and the finite upper bounds of the time integrals, everything is identical to 
the calculation of the propagator of a scalar field ip interacting with a scalar -0 with vertex c/?^^ 
in Minkowski space. In particular, the self-energies S> and S< are related to the time-ordered 
self-energy by 

T.F{x-y) = 0(x°-y°)S>(x,y)+e(y°-xO)S<(x,y) (84) 

and the counter Lagrangian is obtained from the residue of the pole (in dimensional regulariza- 
tion) of the time ordered version of A^rcd 

In order to find the correct combination of curvature tensors we need to identify the tensor 

structure giving the term ^p^p'^^ . This is where the remark made at the begining of this section 

becomes useful. The vertex is nothing but — ^ Jd'^x^/—g 6gabT°-'^ (with /iqo = hij=i = (f>i given 
by (j73p ). Since the tensor structure is the same as in the Lagrangian formalism, the counterterms 
must also be the same. We therefore conclude that the corresponding self-energy operator is 

/ 2\ _ f ^'^P ~^ ~^ 



where 



^l^pua ^ ^pp^ua _ Ij^M^^P^ (87) 

which has already been calculated in e.g. l^, 14, 15, 1^, 17 1. The counterterm is 

Ar = !(fxCabcd{x)C''^"^{x) (88) 

1207r2 (47r)2 4-^7 ^ ) \ ) \ J 

in addition to B? (see the comment (2) in sec. IIIIBp . To recapitulate, this conclusion was 
reached after noticing that the analytic and tensor structure of the loop are identitical in the 
Hamiltonian and Lagrangian formalism. The counterterm must therefore be the same. 

We can check this by a direct calculation. Notice that (j54p depends only on the combination 
ep'C,' as and is therefore adequate to renormalize its contribution. Combining (177p with (155p . 
we get 



(89) 



The dimensionaly regularized integral J is 



which fixes the value of Cw to 
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VI. CORRECTION FROM FOURTH ORDER INTERACTIONS 

I now turn to the corrections from the fourth order interactions 

A4P((/) = -2Im f dh (//4(tl)Cq(T)C-q(r)) (92) 

I first consider the effective fourth order interactions of the previous section (the last two terms 
of ([TOj) ). They turn out to be problematic. I continue with the terms ([50]) . I show that: 1) the 
contribution of (|5ip vanishes. This means that first solving the constraints and then quantizing 
the theory introduces no correction at one loop, which is a non trivial result. 2) So do the 
corrections from (j52p . 3) The remaining terms are renormalized by the same counterterms as 
the tadpole. Contrary to the latter, a secular dependence ln(— gr) appears. 



A. Remaining terms of eq. (j70p 



Let us begin with the term 



2Re J^dn {{TiHs{ti)Q{t)) - {Q{T)TiHs{n))} 

T 

dn p (WgCg[Cg(T)]^ - C.C.) a''(TooCj'^) 



2— Re 



(93) 



When not specified, the argument of the functions under the integral sign is t\. Since the part 
coming from the Wick contractions of the C's is purely imarginary, we only need to calculate 
the imaginary part of (Toocr'^). The result is 



-e(3-26)|C; 



0:4 'hisi) 



with the divergent intergral 



pk, 



J2{q) = Jd^pd^k {2nf6{(i + p + k) ^{p + k) 

The divergent coefficient has the correct dependence to be absorbed into ([60]) . 
Things do not work so well for the other two interactions. The second term is 



(94) 



(95) 



Re J dn {iVg {iv'g + 2p^) [C*(r)]2 - c.c.} a^(rooroo)(q) 

+ Re l^dn {^,2[c*(t)]2-c.c.} a^(roor^o)(q) 



(96) 



Since {TqqTqq) is real, we are left with the second term. Let us write Im{TQoTQQ){q) = Kq + TiK2 
with the (^-dependent divergent integrals 



= j k°lV2l'{2A:V(A:+p) + ^p-k(p.k-/)| 
K2 = f |a°|V^lW{p-k(p-k-/)+/} 



The measure is the same as in (19511. We are left with 



Im j dn u:la\ {Kq + 7^X2) 



(97) 



(98) 
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to calculate. The integrals Kq and K2 have respective dimensions q'^ and g^. After the time 
integration, they are multiplied by q~^, giving thus a correction j4g|(^g|^. The divergent coeffi- 
cient A cannot be renormalized by any of the counterterms. We face the same difficulty with 
the last term 

mt), [J^{t), Qim) ^ Re{u;gC; {u;,C;-u;;Q] X a\T)Uq,T) (99) 
The part depending on C, is q^^ (^ICgP)^5 ^"^^ ttis part depending on matter is 

^3(9, t)= j d^i...cPk2 <5(q + Pi + ki) 5(q + p2 + ks) (roo(pi, ki)roo(p2, k2))(r) (100) 

We stumble against the same nonrenormalizable divergencies as with the tadpole (j65p . for 
which the approximation ()13p was bearing the responsibility. We can reasonnably assume that 
both problems have the same solution. What distinguishes [J-, H3] from the other two is that it 
is the only one with H3. As for (j74p . the derivative couplings make things inside the loop look 
like we were calculating on a Minkowski background, for which we do not need to invoque the 
adiabatic expansion. 



B. No gauge anomaly 

I show that -fT" does not contribute to the power spectrum at one loop. I begin with an 
elementary remark. The terms of the form 

{&d,a) V' , (adia) V' (101) 

have a vanishing contribution in the vacuum because of rotational symmetry. Several terms can 
be cast in this form after an integration by parts over the spatial coordinates: 

diwi , Wagam , -Ad^B d^B^ dtC , 2 ^3C - ^ d.d^B diwl, (102) 

This leaves us with three terms, all containing a Bm- I now show that their contribution also 
vanishes in the vacuum. Consider for instance d^BmdiB diC, which gives 

y* dh a^ih) Jd^p p2 (i?„(p)) Jd^m' 1 • 1' imr, CqC-q]) (103) 

Refering to expressions (j22|) and ([23]) . the Fourier transform of B^ is 



k • k' 

4/3p2(5„,(p)) = 6{k + k' - p)(<jk<Tk' ^-CTkCTk' + Wamip)) 

„2 



jd^k (^Idfcp + ^kfcp) + 4y(a^(k = 0)) 



(104) 



so that the integral fd^p p^ (i?m(p)) is just the term of (|104p inside the brackets. To see that 
this term vanishes identically, turn to the energy constraint at second order 

-6C^ + 12p diB diC + 4C d^B + Ap d'^B2 + didjB didjB - d^B d^B + + 



-8a^ - 8C ^ + 2^^ - 2aV - 4a2l/ (105) 
a'^ 
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and take it's Fourier transform, followed by the limit g — > 0. The result is 

r,2 



4V{a^{q = 0)) = - Jd^p (^\ap\^ + ^kpp) (106) 



which proves the announced result. The contributions of y3C — C/pj d^Bd^Bm evidently van- 
ishes for the same reason, and so does the contribution of (S^ — C,/ p) didjB didjBm because 



jd^p p,pj (BM) = ^-f jd'p P' {BM) (107) 



In conclusion, i/f gives a vanishing contribution at one loop. 



C. Remaining contributions of and trace anomaly 



The terms -ffcxtra do not contribute. To see this, first notice that the expectation value of the 
matter fields produces a (5(1 + 1') which combined with the other Dirac distributions of eq. ()52|) 
yields (5(q). Second, recall from the discussion below equation (j39p that the coefficients C, D, ... 
vanish for the zero mode of C,. This implies that the whole contribution of Hgxtra. vanishes at 
one loop. 

The treatment of the remaining corrections is very similar to the renormalization of the 
tadpole. It pays off to write the remaining terms of H4 exhibiting Tab, 

A4P = -Im I [CirTf j^dn 3Cy{Too) + 2(C,' - C,a,)a^5'^T,j)^ (108) 

After substitution of the expression (j67p of the renormalized energy-momentum tensor and 
integration over time, one finds 

A4Pren = "8/3 (?3|C°|^ { ln(-(?T) + cte + 0(-(?t)} (109) 

We note that this correction is identical to the one of the counterterm p'^ which is the last 
counterterm that had not yet found a use. 



VII. DISCUSSION 
A. Summary of the results 

The technical point treated in this work is the renormalization of the correlation function of 
from matter loops in the Hamiltonian formalism. The reasons for this are explained in sec. [H 
At one loop order contribute two types of terms, the ffi'st with two insertions of the three vertex 
(aa (I49p . the second with one insertion of the four vertex CC'^f (|50|) . 

But with one exception, I checked that the divergences can be renormalized by the same 
counterterms as in the Lagrangian formalism. I argued in section IVI Al that this exception is a 
technical problem, not a fundamental one, and I leave its resolution to future work. I showed that 
the tadpole can be consistently renormalized to zero. Concerning the power spectrum, I showed 
that the corrections from the three vertex are renormalized by the square of the Weyl tensor. 
Concerning the corrections from fourth order interactions, I showed that the terms generated 
by fixing the gauge and inverting the relation between fields and their conjugate momenta have 
a vanishing contribution, and I finally showed that the divergences can be absorbed in the 
counterterms E4, and OR. 
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We can finally stage the final result, and make contact with the introduction. The relative 
correction to the power spectrum at late time is 

^ = -B eGH^ In (j^^ - AGH^ In (-gr) + 0{q^) (110) 

where A and B are constants and fi is the arbitrary renormalization scale. The second term 
comes from an asymptotic expansion, valid for qr <^ 1. The invariance of the action under 
dilatation was identified as the symmetry principle responsible for the conservation law of 
We have now verified that at one loop, the trace anomaly can be the cause of a violation of this 
invariance. 

The following three comments concern the ;S-term. They aim mostly at showing the physical 
transparency of the Hamiltonian formalism. In particular, comments IVIIBI and I VHP I are not 
new but make contact with well known results obtained in the Lagrangian formalism. The 
comment IVII El concerns the choice of the approximation of the modes and its incidence on the 
^-term responsible for the violation. Section IVIIFI is an attempt to open up perspectives on 
the backreaction problem. 



B. Physical interpretation of the counterterms 

Two of the counterterms were found to be 

JV^C^^'^Cabcd , [R'''''''Rabcd-4R'''Rab + R^) (HI) 

in addition to R^ and DR. They are associated with topologically distinct graphs giving ana- 
lytically distinct corrections, eg^ln(g//i) and respectively. This is in perfect agreement with 
what we know about the relation of these counterterms to the trace anomaly. Let me briefly 
recall these facts. 

The square of the Weyl tensor and the Euler density indeed have two different geometrical 



and physical meanings. In the nomenclature of [20|], the former is called the type B anomaly 
and corresponds to an effective action which is not invariant under dilatations since it depends 
on the renormalization scale /i 

^ fd'x^C^'^^lnf^r^Cabcd (112) 



38407r2 V ;^^2 

The Gauss-Bonnet term (type A anomaly) on the other hand is the unique integrated scalar 
density (quadratic in the curvature) which is invariant by dilatation. We understand in those 
terms that renormalizes the correction from (I7ip . since it is the only origin of the ln{q/ fj,), 
while the Gauss-Bonnet term renormalizes the polynomial corrections from the fourth order 
interactions. 



C. Further role of the energy-momentum tensor 

The ;S-type corrections oc (7^1n(g/^) from conformal fields (scalar, Dirac fermions and gauge 
vector fields) have been calculated in 0] and 21]. The stricking result is that they only differ 
from each other and from the case of a minimally coupled scalar by the value of the constant B. 
Moreover, this constant has a universal sign, to wit it is positive for every type of matter field. 
I explain in turn in this section and the following one the physical mechanism behind these two 
results. 
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The universal form of the i3-corrections is really non trivial because of the different time 
behaviours of the mode functions of minimaly and conformaly coupled fields. Conformal modes 
(jl5p redshift, i.e. their amplitude oc a~^, while minimaly coupled fields (jl4p are parametrically 
amplified and asymptote to a constant (this is how the primordial spectra are generated by 
inflation). Naively, one would therefore expect loop corrections from the latter to be enhanced 
by the amplification factor {aH/q)^. 

This is not so because, and only because, the metric perturbations couple to the energy 
momentum tensor. This is revealed by an inspection of Hj,, 

H3 = a* j dT<Pi (TOO + T\) ^ ea^C^e"'^" + T\) (n) (113) 

I used the expression (|13p and I implicitly made the Wick contractions of Q with the external 
leg carrying a wavenumber q. 

Let us first consider conformal fields, and let us take it to be a scalar for the illustration. 
Using (fT2|l . we have 

(rOO + r.)(p,r) = |ap + papp + ^|ap|2 (114) 

where p is the momentum circulating in the loop. Using the solution (|15p of the mode equation, 
we see that the combination 

d^ + Ha^ = i^a^{t) = i^ale-'P^^ (115) 
scales like the gradient term in (|114p . so that 

(T°° + n)(p) = |-^a° (116) 

o Qj 

The factors of a in (jll3p therefore cancel. The only remaining dependence on the background 
is through the slow-roll parameter, and the only dependence on ri and T2 comes from the Wick 
contraction of cjp(ri) from a first vertex with (3"p'(r2) from the second vertex, producing an 
additional term e~^^'P^^ see Eqs. (|77ti80p . Hence the ri and T2 appear only in phases, 

making the time integrals finite for r ^ 0. This is Weinberg's theorem in action (see section |T]), 
and this is what is expected of conformal fields. 

Let us now turn to the case of a minimaly coupled scalar field. Then, 

(TOO + r,)(p,r) = 2|ap|2 (117) 

This quadratic form is very different from (I114p but we see that the solution of the mode equation 
gives, 

<^p(ti) = -^crOe-^^^i (118) 

and comparing with (jlSp and (jllSp . we find that the factor ap/p of the amplitude that 
is gained by the amplification is exactly compensated by an extra factor p/ap from the time 
derivative. with minimaly and conformaly coupled fields are therefore of the same order 

m = 0] ~ m = \] - epc°e"*^" (119) 

They only differ throught their functional dependence on p and p'. 

For this to happen, we see that having derivative couplings is necessary but not sufficient, 
as illustrated by the difference between (jll4p and (jll7p . The coupling of the metric to the 
energy-momentum tensor is instrumental. 
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D. Unitarity 

Prom the expression (|112p of the non local part of the effective action, we immediately obtain 
its imaginary part 

since ln(— = ln(|(/^//x^|) — iiT9{—q'^). Using the argument of [l^ based on the spectral 
representation and that I do not reproduce, we can now explain the fact that the constant B 
in (jllOp is always positive, independently of the spin and the coupling: its positivity means the 
positivity of the spectral density. In other words, ;S > by unitarity. 

E. A comment on the de Sitter approximation 

I now discuss the de Sitter approximation employed in the previous calculations. It is used 
in two instances. One is the renormalized value of the energy momentum tensor (|67p . the other 
is the approximation for the modes (jl3p . Jointly, they are responsible for the secular terms 
ln(— gr) in equation (|lU9p . 

As explained in section ITVl the approximation (jl3p is inappropriate to renormalize the energy- 
momentum tensor. It is necessary for the renormalization of composite operators, which are by 
definition a product of fields in the coincidence limit, to use an adiabatic expansion. By contrast 
(jl3p is only valid around the time of horizon exit and does not capture the correct ultraviolet 
structure of the vacuum. So why is it that we seemingly did not need to use the adiabatic 
expansion to the renormalization the corrections from in section |Vl and can we trust the 
result (j89p ? I think that we can, and I attribute the success of the approximation (jlSp in 
that case to the fact that describes the conformal sector of the theory. I mean this in the 
sense that the counterterm is the conformaly invariant density y/—gC'^, that H-^ is identical for 
both minimaly and conformal fields (see comment IVII Cp . and that contrary to J- and H4, it is 
suppressed by a factor of e which measures the deviation of the background from the conformaly 
invariant de Sitter space. These unique properties of are what made a direct identification 
of the counterterm possible, because inside the loop, everthing looks like in Minkowski space. 

I also point to another limitation of (jlSp . at low wavenumbers this time. I calculated the late 
time expression of (jl09p in a model of power law inflation. The only modification with respect 
to the de Sitter approximation is the expression of the modes, because the renormalized value 
of the energy momentum tensor still scales like l/r^ as in the de Sitter limit (this is because in 
both case the scale factor is a monomial in 1/t). The details are presented in appendix [Cl and 
I quote the result of the late time asymptotic expansion, e.g. 

e \apj 

It is enhanced compared to the classical correction (q/ap)'^ (on the r.h.s. of ([3])), but it becomes 
negligeable after 0{e~^) efolds following the horizon crossing of the mode. So in that case we 
may conclude that the correction is finite. 

The outcome of this superficial analysis casts some doubts on the reality of the secular term 
(jlOgp . which now appears to be an artifact of the approximation (jlSp . The asymptotic decay 
of A4'P in (jl2ip is the consequence of having a red spectrum, i.e. T'{q) ~ g'^s-i with a spectral 
index ns < 1. I have two reasons for favouring this case. One is that it is the scenario singled 
out by observations The other reason is more speculative: because a factor of q always 

appears with a factor 1/a, a spectrum with a blue tilt over a large number of decades produces 
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corrections scaling like some positive power of a{t). In that case we cannot take the limit r — > 
at any stage of our calculations. Although the connection between the power spectrum of C 
and backreaction is not clear (see further comments in section IVIIFp , one could presume that 
this behaviour is a sign of strong backreaction. If one accepts this premise, one of two scenarios 
can happen. One is that this backreaction tends to increase the blue tilt even more. In that 
unfortunate case the game is over until we can solve the full nonlinear problem. The other 
alternative is that the backreaction tends to soften the blue tilt, perhaps enough to turn it into 
a red one. In this happy turn of event we can use in first approximation a red spectrum over 
the whole range of q. Doing so, we certainly make a quantitative error on the amplitude of the 
correction, but not a qualitative one. 



F. Perspectives on the backreaction problem 

This brings me to my last comment. As we learn to better handel the calculations in the 
Hamiltonian formalism, we can contemplate tackling the backreaction problem with these tools. 
One of the foremost advantages of the formalism is its physical transparency, since we work 
directly with the physical degrees of freedom and the two helicities of the graviton) and use 
a physical clock (the inflaton field is equal to its classical value). If need was, I hope to have 
convinced the reader of this with the first three comments. 

The theoretical obstacle to the backreaction problem is the identification of the relevant 
observables to quantify the effect. Most of the efforts have been invested in the calculation of 
the impact of non linearities on some measure of the local expansion rate, with bewildering 
difficulties to desentangle the physical effect from a gauge artefact ^. 

The expansion rate is but on facet of gravity, and there is bound to be something to be learned 
from observables sensitive more specifically to other effects, for instance to the decelaration 
parameter or to tidal forces. I invite the reader interested in these issues to read the lucid 
discussion of Tsamis and Woodard on this point [2^. I have little to add to this, but for the 
following two remarks. First, the tilt (including the running) of the power spectrum of scalar 
and tensor perturbations is an interesting quantity. A tilt indeed measures the variation of the 
amplitude of the modes at horizon exit in response to a variation of the Hubble rate, via the slow 
roll parameters. In other words, it measures an acceleration. The scale dependent corrections 
in ()110p can indeed be interpreted as a correction to the spectral index 

An, = = -BeGH' (122) 

fling 

and it is interesting that unitary, i.e. B > 0, implies that the backreaction from matter fields 
renders the spectrum more red ^. 



^ If we limit ourselves to the topic ofthe backreaction during inflation, a selected list of references starts with the 
works of Tsamis and Woodard and of Mukhanov & al. [2^ , respectively on the backreaction of gravitons 
and of superhorizon scalar perturbations. Their methodology was seriously criticized by Unruh [2^. Subsequent 
works taking these attacks into careful consideration are those of Abramo and Woodard Geshnizjani and 
Brandenberger [2^, and Finelli & al. [2^. Although these groups consider different observables, their results 
provide evidence that, as Unruh anticipated, the backreaction of fluctuations on the homogeneous mode cannot 
be observed by local observers who should interpret superhorizon adiabatic perturbations as a scaling of the 
coordinates, provided that the adiabatic mode decouples, i.e. provided that the non adiabatic pressure is 
negligible, so that the inflaton is the only physical clock at their disposal. This is also the intuition at the 
basis of ([iJS]). It implies concurrently that in single field inflation, both backreaction and a violation of the 
conservation law of can only be sourced by derivative couplings, in agreement with eq. Q. 
The corrections to the slow roll parameters from the fluctuations of the inflaton have been calculated by three 
different groups [sil. Is^. l33l| . They each make different simplifying assumptions, employ different methods of 
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The second observation is that it seems reasonably feasable to calculate a quantity such as 
the expectation value of the squared Weyl tensor. The latter is a sum of n-point functions of C 
and gravitons 7, 

{C'^'^'Cabcd) = 32a-\epf{C'iT, x)C'(r, y)) + Oi^f^X') (123) 

This observable has direct physical significance since it is a measure of the curvature. Moreover 
it is a scalar which vanishes on a flat Fiedmann Robertson Walker background, meaning that it 
is a gauge invariant measure of the deviation from this background ^. 

There are several technical difficulties that need to be resolved before undertaking this cal- 
culation. We have already been made aware of one of them, namely the approximation used for 
the modes, both in the infrared and ultraviolet. This is nothing new however. If we go as far 



as two loops, as is necessary in de Sitter 2j|, we will have to learn how to handle overlapping 
divergences in the Hamiltonian formalism. 
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APPENDIX A: CHANGE TO THE LONGITUDINAL GAUGE 



All the quantities in this section are first order. Scalar pertubations are therefore decoupled 
from vector and tensor modes, in any gauge. Consider the following parametrization of the 
metric and inflaton 

ds"^ = a^{T) {- (1 + 20) dr^ + 2diLodTdx' + [(1 - 2^^) 6ij + 2didjE] dx'dx^} 

ip{t,ji.) =if{t) + 5ip{t,^) (Al) 

In a first order coordinate transformation parametrized by 

T^f = T + f, x'^i' = x' + diC (A2) 

the five scalar perturbations defined above vary as 

a a 
6^ = 6ip- ^{t)f (A3) 

By definition, in the comoving and longitudinal gauges we have respectively 

E, = 0, 5ip,{t,^) = 

Ei = 0, Bi = 0. (A4) 



calculation, and obtain different results. In particular, one source of discrepancy can be traced back to their 
use of different gauges. Contrary to (|6]), these gauges do not use the inflaton as a clock, which makes difficult 
to assess the reality of the effect to local observers. 

Note that the contribution of the scalar perturbations is 0{e) (because of the normalization of the modes p3p 
contains a factor e~^), so that the leading contribution comes from gravitons. This contrasts with the intriguing 
results of Losic and Unruh [s^l, to wit, certain gauge invariant combinations of the scalar perturbations are 
dominated by second order fluctuations. Clarifying such discrepancies should prove to be fruitful. 
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The passage from the comoving gauge to the longitudinal gauge is therefore realized by the 
coordinate transformation 



i.e. = —uJc and ^ = 0. With (j)c = ipc = — Ci = —j/ + ed '^C,' , and the linear equation 



Tc ^ n = Tc- UJc, X\= xl (A5) 

= —i^^ and ^ — n WifV, a, — ,h — — j- e/)-2/^' 

([TT]1 . one obtains 

(Pi = 1^1 = -ep'd-\' = -epa'd-^C (A6) 

The l.h.s. of this equality is recognized as the term multiplying 20"'^ = d"'^Tab (for a minimaly 
coupled field) in (j7ip . In words, the vertex is the interaction between a and the metric in the 
longitudinal gauge. This can also be shown explicitely. In a coordinate transform x'^ i— > + ^'^, 
the action varies by 6^S = — jd'^x^J—g ^°''VhT\, which induces the following correction of the 
C(7C7-vertex, ^ 



Jd^x^ I ^T°° - 2diUjT°' - 2CT\ - lo {dtT^^ + SpT"" + diT^' + pT\) I (A7) 



Integrating by parts the second, fourth and fifth term and using the linear equation (|lip . one 
finds 

Si[Caa] = Jd^x^cPi (TOO ^ ^^.^ _ J^^^g^ (a^T^'^Lo) - jd\di {a^T^'u;) 

= - jdt (^H3 + dtJ"^ ~ j'^^ j (a^r°'w) (A8) 

The last term is a spatial gradient and therefore does not contribute. The term T is now 
clearly identified with the generator of the (first order) canonical transformation relating the 
longitudinal gauge to ([7]). Indeed, the Einstein Hilbert Lagrangian density changes by a total 
derivative in a gauge transformation. So does the Lagrangian of any system in a canonical 
transformation. So a gauge transformation can be seen as a canonical transformation. 



APPENDIX B: EXPRESSIONS OF THE GEOMETRIC COUNTERTERMS 

1. Weyl tensor 

The first counterterm is the square of the Weyl tensor. Since the latter vanishes for a 
flat Priedmann Robertson Walker space-time, we only need to calculate it at first order. The 
components of the Riemann tensor are 

R\ = {P"(l + 20 + ep'C] + d,d, {ep'd-\') (Bla) 

R%k = 0{e) (Bib) 
^'ifc, = (p')'(l + 2C) {6\h-5'^5,u) 

+ [Sjididk + Sikdjdi - dikdidj - Sijdkdi] (ep'5~^C') (Blc) 



® The passage to the longitudinal gauge generates of course other terms that I have not discussed here: corrections 
to the pure gravity vertices, as well as interactions with the inflaton perturbations. 
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from which we obtain the components of the Ricci tensor and the Ricci scalar respectively 

^00 = -3p" + (2e - 6)p'C' (B2a) 
Ro^ = 1Zd,uj = [p" + 2{p'f] diUJ (B2b) 
Rij = 1Z{l + 2Q)5ij (B2c) 
a" 

a^R = 6— + Aep'C (B2d) 
a 

The components of the Weyl tensor in four dimensions are 

C%cd = R%cd - \ Rhd - Rbc - Qbc R^d + 9M R^c) + I 9bd - 6^d 9bc) (B3) 
and I find 

C%^ = ep'D,,{d-\') (B4a) 

C%k = 0{e) (B4b) 

C'^ki = [5,1 Afc + Dji - 5,k Du - 5u D,k] {ep'd~^C') (B4c) 
I introduced the notation 

Dij = didj - 5,jd^ (B5) 

One can verify that these components have the correct symmetry and that the identity C^q^ = 
—C\kj implied by the vanishing of the trace is satisfied. Note that eqs. (|B4p depend only on 
the quantity (see appendix |A|) 

= -ep'd-'^C (B6) 

and have therefore the desired form to serve as counterterm for the vertex H^. Finally, the 
squared Weyl tensor is 

V^gC' = Ne'^{A{C%^)'-A{C%,f+{C),,)'} 

= 8{ep'f [didjd-\' d,d,d-\' + 3(C')'] + O(C') (B7) 
or in the momentum representation 



(B8) 



with Cw a constant. 



2. Second order expressions of the Ricci tensor and Ricci scalar 

For the other counterterms we need the expressions of the Ricci tensor at second order. I give 
for reference the expressions of the connection coefficients F^^ = \g°''^ [dbQcd + dcQbd — ddgbc\'- 

FOq = + a' + - aa' + diadico + p'{dujf (B9a) 
Fqq = di [a + uj' + p'uj + 02 + "^2 + + p'wi 

+ (2C' - a')diUJ + {a- 2C)dia - didjiodjcu (B9b) 
Fqj = di [a + uj' + a2 + ^2) + p'wi — ('diUJ — adiU — didjUjdjUj (B9c) 
r?,- = {p' - C + 2p'C - 2p'a2 + 2p'e + p'a^ - 2CC' - diCd.oj] 

- didj {uj + 0J2) {diWj + djWi) + 2 (a - C) didju (B9d) 

Toj = {p' + C) ^ij - diUjdja - p'diUjdjUJ + ^ {djWi - diWj) (B9e) 
Ff^. = diC 5) + d,C 5\ - 6i, {dk (C + p'^ + P'w2) + p'wk + (2p'C - OdkC] (B9f) 
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The second order expression of the Ricci tensor are: 

62R00 = 3p'a2 + 5^ ["2 + 4 + - {diof] + nidujf - SiC'f 
+ SdiC'ditu + Sp'diCdiUj + SdiCdiu;' + diCdiU 
+ (2C' - a') + (a - 2C) S^a 

I separate Rij into diagonal and non-diagonal parts according to 

Rij = 6ij + Rf^ 



with respectively 

= e'^<n [1 + q2 - 202] - p'a'2 - p'd'^LJ2 - 3(C')^ - (.dCf 
+ 2ep'CC - diCd.a + (C' - 2p'C) d^uj - 5p'diCd,uj - [diCdM' 

Rf^ = - didj [02 + uj'2 + y^'A - \ {^i^'o + ^i^O - i^^^j + 

+ adidja + (C' + a' - ip'C) didjUJ + 2 (q - C) didjco' 
+ diCdjO + djCdia + 5iC9jC 

We do not need the second order expression of Roi because it is a first order quantity which 
multiplies other first order quantities: to calculate the Ricci scalar g"'^Rab, it is multiplied by 
g^^ oc diLO, and to calculate R°'^Rab it is multiplied by i?^*. The second order part of the Ricci 
scalar is 



(BIO) 
(Bll) 

(B12a) 
(B12b) 



a'^52R = -12—02 - Qp'a'2 - 2(9^ 
a 



a2+uj2 + 3p'u!2 - -jiduj)"^ 



- 6 [diCdiUj^ - ISp'diCdiLu - 2diCdia 

+ 2Cd^a + 2{C' + a' - bp'C) d'^uj + 2 (a - C) 9^ 



8e(C')' + 6ep'CC'-2(aC)2 
(B13) 



3. The remaining counterterms 

Using the solutions of the constraints and the linear mode equations to simplify quadratic 
terms containing a C," or cr", and dropping spatial divergences, I finally get 



2 

Cr2 



(B14) 



1 



R^ {l + 3C + a + lc^ + 3Ca) + SRep'C + R ( {2 + e)Too + ^ 



+ iR -^{d'CY - "5 a + (1 + 5e + eg) ad\ - (3 + 2e) C^^C 



+ 36p'\c' + 624e(p'C')^ - 296(ep'C') 
I used the shorthand notation R = 6—, and 2Tqq = (cr')^ + (9cr)^. I neglected terms of order 
higher than two in the slow roll parameters e = 1 — and £2 = 

I give the expressions of the other counterterms only at first order in C,. Several combinations 
are interesting. The integrated Euler density reads 

Jd'^xV^E^ = jd^X^ {R'^'^'^Rabcd - ^.R'^'Rab + R"} 

= 8 Jd^xNe^< {3pVWeC' + 0(C')} 

= 8 Jd^x{-dr{p'')+p'^C' + 0{e)} (B15) 



„V\2 
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and I used p" = /o'^(l — e) which fohows from the Friedmann equations. The hnear term 
coming from the measure has been ehminated by an integration by parts. The remaing term is 
proportional to eC' and vanishes in the de Sitter limit e — > 0. One can verify that at linear order 
the combination 

\^9{r\R\-\r^ = 4 Jd'x{dAp")-p"eC' + 0{e)} (B16) 
is in agreement with the identity 

E4 = - 2 (^R\R\ - ^R^^ (B17) 

which follows from (|B20p below. The last geometrical counterterm which is fourth order in the 
derivatives is f\/—g OR- It is a bit ambiguous since it is a total derivative. If we drop the spatial 
divergences alltogether and expand the remainder at linear order, we obtain 

'^OR = Jd^xda[V^g''%R 

d^xdr{^R! 

= -24 jd^xe^p'\' (B18) 

It is subleading in the first slow roll parameter. If we do not drop the spatial derivatives, two 
additional linear terms fd'^x e^a^^d^R + R'd'^oo) are included. We can integrate them by parts 
to produce a second order term oc diQdiC,' . 

Finally, it is easy to show that the counterterm ()B14|) is comprised of all the terms at most 
quadratic in C, with a total of four derivatives and a maximum number of two time derivatives 
of Q. Consider first 

^b^a — \^0) + qU i + U jK ■ 

= {R\f + ?,'R? + 2lZ5'^52R,j+0{C^) (B19) 

The product of the 0-z terms is third order because i?*Q oc diu) and R^^ = 2TZ{a + C,)diUJ, and 
the linear term of i?*^ = TZ5ij + 52Rij vanishes. We see that at second order ()B19P differs from 

R^ = (yR\ + 'i'R- + 5^^52Rij)'^ only by a different combination of the same terms. Similarly for 

R''''"^Rabcd = C^ + 2R'''Rab - ^R^ (B20) 

which follows from ()B3p . The last counterterm J \J—g Di? does not contain any new quadratic 
term either, since every time that a C" or a" appears we can substitute the linear field equations. 

APPENDIX C: TIME INTEGRALS FOR A RED SPECTRUM 



a(r) = f— , dla>^ ^ a>0 (CI) 



Consider the model of power law inflation, 

r ' 

The Hubble rate is -^(t) = {^^^ with = ^-^id the only nonvanishing slow roll 
parameter is e = The solutions of the mode equations 
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have the same behaviour as (|13til4p for gr— > -oo. Ti.\P is the Hankel function of the first kind 
of index 

v = l + a>^- (C3) 

I wiU assume that a does not take the values 5 or | for which the discussion below must 
be amended. These special values correspond to e = 1/3 and e = 5/7 which are outside 
the slow-roll regime anyway. To be precise, say a < ^ in the following. Using the identity 
^ = Hu-iiz) - fW^(z), we have 

4 = -f ^^W«(-.r) (C4) 

and a similar expression for dp. The normalization constants and AT/ are fixed by the 
canonical commutation relations, e.g. 



[Cq(t), 7r_q(t)] =i = ^ i^CgC; - C.C 



, - ^ - C.C. (C5) 

\ / z=—gr 

The term in parenthesis is the Wronskian of the Hankel functions and is equal to —4i/{7rz). We 
deduce 

l<l' = J^, \M;;f = l (C6) 

To analyze the late time behaviour of the various integrals, I use the asymptotic expansion 
of the Hankels for small arguments 



n,{z) 



sin(7ri/) 



r(i-i/) \z 



which gives the following expansion of the mode 

C,(r)=C°(l + m,(gr)3+2- + ...) (C8) 

where = \au\ exp(i7ra). The power spectrum is 

q%{r)\' ^ q-'"- (C9) 

Given the condition (jCip . this is a red spectrum (more power in for large wavelengths). One 
can explicitly check that this does not change the contribution of the vertex of section |Vl 

1^2 / /\l-2a 

V.V.ZI^M;, ^ ^(-ro)— [(-.r)3+- + (CIO) 

The limit a ^ continuously converges to (j83p and ()80p . The integrals over time of (jClOP 
converge for r — > 0. By continuity, we can conclude that for small values of a, AV( is given by 
its value (jllOp plus a term oc ar" with a' > that goes to zero for r — > 0. 

To evaluate the integrals ([58]) . ([63]) . (fMI) and (I108p . we need too more ingredients. One is 
the asymptotic expansion of (' / p' 

= -(1 + a)q'r' C° (l + ^b.{qT)'+''' + ...) (Cll) 
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with the coefficient 

6^ = |6^|e^™, = (3 + 2a)(l + 2a) (C12) 

\au\ 

We also need the expression of the renormahzed energy momentum tensor. It can be found in 
[l9| as a polynome of p' and its derivatives. I only quote the result 

"'^'^''^ = ^ (119 + 482a + 0(a2)) 

1 1 / 4S3 \ 

"'^^^^^ = "960^ + — " + ^(«'V (^'') 

Because the scale factor is a monomial of l/r in power law and de Sitter inflation, p' oc 1/r in 
both cases. This means that the renormahzed energy- momentum tensor contributes the same 
to the time integrals in both cases. 

I now use these expressions to calculate the value of the tadpole ()64p . The terms 0{a^) 
of (jClSP give back ([68]) . and the terms 0{a) are multiplied by a factor produced by the 
integration over time. The leading term of (Cq(''")) is therefore 0{a^). The final result is 

(Cq(T)) = #<5(q) X q^qrr (C14) 

where 7^ stands for a silly number. We have exactly the same result for (j63p . 

I then apply (ICSP to the integral (f58]) . and (jl09|) . which are also subject to identical changes. 
For instance 

Im fdnp"{C{n)C'{n)e{T)) ~ q^l' x ^ (^)'" (C15) 

The r.h.s. is approximated by + ln(— gr) for time intervals qr <^ a^^ ~ and asymptotes 
to zero. The large value for ^ 1 is not significant since (jClSP is an asymptotic expansion. 



APPENDIX D: THE CTP FORMALIM 

For time dependent problems, where there is no stable ground state, there exists a covariant 
formalism called the Closed Time Path formalism adapted to calculate expectation values such 
as Green functions. I briefly introduce this formalism for the readers convenience and refer to 



3j] for other short reviews and further references. The important point for us is the structure 



of the one loop correction to the two point function given in Sec. ID 31 



1. Motivation 

Let us consider the time-ordered two-point function, defined in the Heisenberg picture by 

G{x, y) = {^Q\T^H{x)^H{y)\'^G) , (Dl) 

where |^o) is the exact (non-degenerate) ground state at an 'initial' time to. Inside the horizon, 
the evolution of the modes is adiabatic and a vacuum state is well defined. The initial time can 
be taken — > —00 and the initial state in noted |^'o) = |Oin). 

In the interaction picture, for the particular time ordering > ip, Eq. (|Dip becomes 

G{x,y) = (0in|C/7(-oo,x°)(^/(x)C//(x0,y0)(^/(y)C//(y0,-oo)|0in). (D2) 
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where Ui is given by 

Ui{t,to) =Texp(^-i J^dt' Hj{t')^ (D3) 

and Hj[ipi,TTi;t] is the interaction Hamiltonian written in terms of the field variables v7/,7r/ in 
the interaction picture. Let us insert a resolution of the identity at the time iout = +00 in the 
last line of ([D2]) 

00 

x°>y°, G{x,y) = ^ (Oin|;7/(-oo, +oo)|nout) 

n=0 

X {nont\Ui{+^,x'>)y,j{x)Ui{x^y'^)^i{y)Ui{y'>,-^)\Om) (D4) 

When the ground state is stable, the matrix element (0in|C/7(— 00, +cxo)|Oout) is a phase 
and the matrix elements (0in|{7/(— 00, +cxo)|n 7^ Oout) vanish identically. One thus recovers 
the familiar expression for the time ordered propagator in the interaction picture G{x, y) = 
(0 in|TC//(oo, +00) (/J7(x)v37(y)|0 in)/(0 in|C//(+oo, — 00)10 in). This simple expression no longer 
holds when the initial state is not stable since then all the terms in the sum contribute. In that 
case, the pertubative expansion of the Green function is obtained directly from (ID2p after ex- 
panding each of the three evolution operators in (jD3|) . The perturbative theory can be described 
by a diagramatic formalism akin to the Feynman rules for scattering matrix elements. 



2. The 'in-in' generating functional and Feynman rules 

The so-called 'in-in' generating functional of connected Green's functions is given by the 
transition amplitude of two 'in'-vacua in the presence of external sources J+ and J_, 



e'^^-^^'-^ ] = j_(0in|0in) 



n=0 

= Tr|Te*^*°- '^*'^+('^)^«('=VmTe~^^-°- (D5) 

where T is the reversed-time ordered product and J+ and J_ are the two classical sources asso- 
ciated with the two branches of evolution, forward and backward in time respectively. Because 
of these two reverse branches, the formalism is called Close Time Path. In the second line the 
resolution of the indentity at an arbitrary time tout was inserted. In the last line, the generating 
functional was written in a form that allows arbitrary initial states (not only the vacuum). Note 
that the operation of taking the trace couples the forward and backward time evolutions. As a 
result, VF[J+, J-] generates four types of connected two-point functions: the usual time-ordered 
propagator 

5W 

G++{x,y) = i{Tip{x)ip{y)) = J j+=j_=o (D6) 

dJ+{x)dJ+[y) 

the reverse times ordered propagator 

SW 
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and the two on-shell two-point functions (the Wightman functions) 

5W 

5W 

G^.{x,y) ^ ^{^{yMx)) = JJ^^^^^^^W=J-=o (D9) 

The latters are not time ordered because they are build with operators coming from different 
branches. 

The path integral representation of the 'in-in' generating function is obtained from (|D5p by 
doubling the fields, one (/?_)_ for the forward branch and one (/?_ for the backward branch. The 



obtention of Feynman rules was described in details in [34i | and [a]. They are similar to the 
'in-out' Feynamn rules, but for a few differences: diagrams are build from the four propragators 
()D6P - (|D9p and two types of vertices of polarity + (the same as the 'in-out' vertices) or — (with 
a relative minus sign). At each vertex, in addition to the conservation of momentum, one can 
attach only the propagators with end lines of the same polarity. For instance, we can have the 

sequence ...G |_(vertex+)G_| (vertex— )G ... 

All the algebraic theorems of the in-out formalism are naturally extended to the in-in for- 
malism. In particular, the double Legendre transform of W^[J+, J_] gives the in-in generating 
functional r[(^+, (pJ\ of 1-particle irreducible vertex functions. The second variational derivative 
of r[(^+, (pJ\ therefore gives the matrix of self-energies: 

r[^+, = 52[^+] - S2W-] -\fd'x jd'y ( ^+ c^- ) ( ) ( ) (DIO) 

where ^2 is the quadratic part of the classical action, e.g. (|28p . The self-energy matrix elements 



are 



^-+{x,y) = ^>{x,y) 

^+-{x,y) = ^<{x,y) 

S++(x,2/) = ^(x°-y°)S>(x,y) + %0-xO)S>(x,2/) 

S„„(x, y) = ^(xO - y°)S<(x, y) + 0(yO - xO)S>(x, y) (Dll) 

The in-in effective action is renormalized by the same counterterms as the in-out effective action, 

because the on-shell self-energy (as well as higher order vertices) is finite and S is the complex 

conjugate of time order self-energy 



3. Generic structure of the one-loop correction to the propagator from H3 



Consider for simplicity a theory with trilinear couplings iptp'^ between two scalar fields. Using 
the in-in Feynman rules, the one-loop correction to the Feynman propagator is 



AGFit,t',^-y) 



d^z, I d'z2 Gl,{x - zi) \-iT.Uzi - Z2)\ Gl4z2 - y) (D12) 



0, and 



where the indices are summed with the metric C++ = +1, c = —1, and = c_.^ 

the matrix elements Gab are given at (|D6h - (ID9p . The self-energy matrix is 

^U^) = Gt{z)Gl{z) (D13) 
The matrix proguct gives the following integrand 

G++(t,ti) {S++Gii(t2,t) - S+_G-+(t2,t)} - {G+-(t,ti)S_+ - G+-(t,ti)S__} G-+(t2,t) 
= G>(t,ti){S++G<(f2,t)-G'>(t2,t)S<}-G<(t,ti){S>G<(t2,t)-S_„G>(t2,t)} 
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To get the second line, I expanded the time ordered propagators, and anticipating on the inte- 
gration over time I dropped the terms proportional to 6{ti — t) and 6{t2 — t). Subtituting these 
expressions into ()D14p . one finally obtains 
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AGi.(t,t,x-y) = Re / dz" / dz', (G>(t', z^) - G<(t', 2°)) 




X [G>(t,^0)S>(zi,z2) - G<(t,z?)S<(zi,z2)] 



(D14) 



(1994). 



